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We present results from a study of X(3872) — >■ -niij/^l) decays produced via exclusive B — >■ 
A'X(3872) decays. We determine the mass to be Mx{3872) = (3871. 84±0.27(stat) ±0.19 (syst)) MeV, 
a 90% CL upper limit on the natural width of Fjf (3872) < 1-2 MeV, the product branching fraction 
B{B+ -^ K+X{3872)) x B(X(3872) -^ tt+tt" J/V') = (8.61 ± 0.82 (stat) ± 0.52 (syst)) x lO-*^, and 
a ratio of branching fractions B{B° -^ K° X (3872)) /B{B+ -^ K+X{3S72)) = 0.50 ± 0.14 (stat) ± 
0.04 (syst). The difference in mass between the X(3872) — >■ tt'^tv' J/ip signals in B'^ and i3° decays 
is AMx(3872) = (-0.69 ± 0.97(stat) ± 0.19(syst)) MeV. A search for a charged partner of the 
X(3872) in the decays B° ^ K' X+ or B+ -^ K°X+, X+ -^ tt+tt°J/iIj resuhed in upper limits 
on the product branching fractions for these processes that are well below expectations for the case 
that the X(3872) is the neutral member of an isospin triplet. In addition, we examine possible 
J^ quantum number assignments for the X(3872) based on comparisons of angular correlations 
between final state particles in X(3872) — > 7r^7r~J/i/> decays with simulated data for J^ values 
of 1''"'' and 2 '". We examine the influence of p-uj interference in the M{-k^-k~) spectrum. The 
analysis is based on a 711 fb~^ data sample that contains 772 million BB meson pairs collected at 
the T(4S') resonance in the Belle detector at the KEKB e^e^ collider. 

PACS numbers: 14.40. Pq, 12.39.Mk, 13.20.He 



INTRODUCTION 

The X(3872) was first observed by Belle as a narrow 
peak in the TT'^Tr~J/^ invariant mass distribution in ex- 
clusive 5+ -^ i^'+TT+TT^,//'!/) decays [l|, l^- It was sub- 
sequently seen in ^/s = 1.96 TeV pp annihilations by 
CDF [3| and DO [4| and its production in B decays was 
confirmed by BaBar [5] . A recent summary of the mea- 
sured properties of the X(3872) is provided in Tables 10 
through 13 of Ref. g. 

The close proximity of the PDG world-average of 



X(3872) mass measurements, 
0.22 MeV 0, to the moo + 



M; 



avg 



I'D' 



= 3871.56 ± 
mass threshold 



(3871.8 ± 0.3 MeV [3]) has engendered speculation that 



the X(3872) might be a loosely bound D^-D 



*o 



olecu- 



/ar state [8[. Theoretical studies of deuteron-like D'^D*^ 
interactions were reported by Tornqvist in 1994, and 
he predicted bound states for J^'~^ values of ^ and 
l"^"*" '95. There has been considerable theoretical interest 
in the X(3872) line shape in its D^D*° decay mode [id|. 
These discussions are constrained by the current un- 
certainty in the natural width of the X(3872) in the 
TT+TT^ J/V' decay channel, which is rj5s:(3872) < 2.3 MeV 
(at the 90% confidence level) [l| . A measurement of the 
natural width in this mode, or an improvement in the 
upper limit on its value, would be useful input to these 



line-shape studies. 

A close correspondence of the tt'^tt^ invariant mass 
distribution to expectations for p — )■ tt+tt^ decays was 
reported by Belle 11 1 and CDF [12|. This, together with 
the observation of the X(3872) — » jJ/ip decay mode by 
both Belle [13| and BaBar [IJ], establishes the charge 
parity of the X(3872) as C = +1. A comprehensive study 
of possible J^^ quantum numbers for the AT (3872) us- 
ing a large sample of X(3872) — > tt+tt^ J/-0 decays was 
performed by CDF [l5|, ll6|; they concluded that only 
the 1++ and 2 '" hypotheses are consistent with data 
and other assignments are ruled out at the 3cr level or 
above. The A (3872) — >• jJ/tp decay process would be 
an allowed El transition for a 1+"'" assignment and a 
suppressed higher multipole for 2 ^; the observation by 
BaBar and Belle of this process favors l"'"^ [17|. However, 
a recent BaBar analysis of the A(3872) — >■ Tr+Tr'Tr" J/?/j 
decay mode showed some preference for a 2 ^ assign- 
ment J18j. Since bound molecular states are predicted 
for J^'-^ = 1++ but not for 2 ^, an unambiguous exper- 
imental determination of the spin-parity of the A (3872) 
is an important input to the understanding of this state. 

Another proposed interpretation for the A (3872) is 
that it is a tightly bound diquark-diantiquark four-quark 
state [19|, in which case two neutral A (3872) states - or- 
thogonal mixtures of cucu and cdcd - are expected to ex- 
ist shifted in mass by 8± 3 MeV. The authors of Ref. [13] 
suggested that these two different states might result in 
different A(3872) masses in the B^ -^ K~^ tt~^ Tr~ J / ip and 
_B" — >■ Ar°7r+7r^ J/V' decay chains. BaBar measured the 
A(3872) properties separately for these two channels and 
found a mass difference (AM = 2.7±1.6±0.4 MeV) that 
is consistent both with zero and the lower range of the 
theoretical prediction [20|. CDF used a comparison of 
their measured A(3872) — >■ TT^Tr~J/ip line width with 
their experimental resolution to establish a 95% CL up- 
per limit of AM < 3.6 MeV, for equal production of 
the two states |2l| . These results are not definitive tests 
of the prediction of Ref. [19| ; the statistical significance 
of the BaBar signal for B° -^ K°X{3872) is marginal 
(9.4±5.2 events) and the interpretation of the CDF limit 
depends upon the unknown relative production strengths 
for the two different states. Thus, a more precise com- 
parison of the A (3872) produced in B+ and -B" decays 
is needed. 

In the diquark-diantiquark scheme, the A (3872) is ex- 
pected to be the /3 = member of an isospin triplet. 
Since the dominant weak interaction process responsi- 
ble for B — > KX{3872) decays is the isospin conserving 
b — >■ CCS transition, the charged I^ — ±1 partner states 
(that decay via A+ — > p'^J/^p) are expected to be pro- 
duced in B decays at a rate that is twice that for the 
neutral A (3872) [22] ■ The BaBar group studied the pro- 
cess B — > KiT^-K^J/tp and placed upper limits on the 
product branching fractions for X+ -^ tt^tt^J/'ijj that 
are below isospin expectations [23|. 



Here we report on a study of A (3872) -^ tt+tt" J/V' de- 
cays produced via the exclusive decay B — >■ X A (3872). 
We use a 711 fb~^ data sample that contains 772 million 
BB pairs collected in the Belle detector at the KEKB 
energy-asymmetric e+e^ collider [2j|. The data were ac- 
cumulated at a center-of-mass system (cms) energy of 
^/s — 10.58 GeV, at the peak of the T(4S') resonance. 
KEKB is described in detail in Ref. (li. 



DETECTOR DESCRIPTION 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer cylindrical drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux- return located outside of the coil is instrumented 
to detect Kl mesons and to identify muons (KLM). The 
detector is described in detail elsewhere 12611. 



B -^ Kii+-K-J/ip EVENT SELECTION 

We select events that contain a J/ip — 7> €+^~ (£+£" = 
e+e~ or ijl^ijl~), either a charged or neutral kaon, and a 
TT+TT^ pair using criteria described in Refs [l[ and [27|. 
The leptons from the J /ip — >■ l^i~ decay are required 
to pass minimal lepton identification criteria and the 
invariant mass of the pair is required to be in the 
ranges -21 MeV < {M^+^,- - mj/,p) < 20 MeV and 
-24 MeV < (Afe+e- - mj/^) < 20 MeV, where ruj/^ = 
3096.92±0.01 MeV is the world-average value for the J ftp 
mass [7|. For J/^p -^ e+e~ candidates, photons within 
50 mrad of the e"^ and/or e~ tracks are included in the 
invariant mass calculation. The number of events with 
multiple J/tp candidates is negligibly small. Candidate 
_ftr+ mesons are charged tracks with a kaon identifica- 
tion likelihood that is higher than that for a pion or a 
proton; neutral kaons are detected in the Ks — > 7r"'"7r~ 
decay channel using the Ks selection criteria described 
in Ref. [28|. The charged pions are required to have a 
pion likelihood greater than that of a kaon or a proton. 
Some events have more than one acceptable combination 
of hadron tracks. In these cases, which include 3% of 
the events in the signal region, the tracks with the best 
vertex fits are used. To reduce the level of e~^e~ — > qq 
(q — u,d, s or c-quark) continuum events in the sample, 
we also require i?2 < 0.4, where i?2 is the normalized 
Fox- Wolfram moment [29|. 

Events that originate from B — > K-K^n^J/ip de- 
cays are identified by the cms energy difference /S.E = 
jTjcms _ E^l^ and the beam-energy-constrained mass 



^bc = v/(^bSrm)' - (Pb")'> where i?^,^^ is the cms 
beam energy, and E"^^ and p'j^^ are the cms energy and 
momentum of the Ktt^tt^J/iP combination. We select 
events with Mbc > 5.20 GeV and -0.15 GeV < AE < 
0.2 GeV. We define signal regions as 5.272 GeV < Mbc < 
5.286 GeV and -0.35 GeV < A^ < 0.03 GeV; these cor- 
respond to ~ ±2.5(7 windows around the central values 
for each variable. 

In addition to selecting B — > KX{3872) events, these 
selection criteria isolate a rather pure sample of i? — >■ 
Kip' , ijj' — >■ iT~^'K~ J/ij} events [30|. These events are used 
as a calibration reaction to determine the Mbc, AE and 
M {tt^ TT~ J / tp) peak positions and resolution values, and 
to validate the Monte Carlo-determined acceptance cal- 
culations. 

For each event we compute Af (tt+tt^ J/^) from the 
relation 



Min+n-J/yj) - M;T:%+^- 



'■J/ipj 



(1) 



where M"°^l,,, 



and M^"^! are the measured 



7r+7r~£+£~ and £'^£~ invariant masses, respectively. For 
studies of the ip' — >■ tt+tt" J/?/" control sample we use 
events in the interval 3.635 GeV < M(7r+7r~ J/f/)) < 
3.735 GeV; for X(3872) studies we use 3.77 GeV < 
M{tt+tt-J/iI)) < 3.97 GeV. The M{TT+n- J/iP) sig- 
nal regions are defined as \M{tt~^tt^J/iP) — AfpoakI < 
0.009 GcV, where Mpcak == 3.686 GeV and 3.872 GeV 
for the tjj' and X(3872), respectively. We select events 
with a dipion invariant mass requirement of M^+^- > 
{M{-K^'K~ J/ip) — {mj/^ -t- 150 MeV)), which corresponds 
to M^+^- > 625 MeV for the X(3872) and > 439 MeV 
for the iJj' events. After this requirement, which re- 
sults in a 6% signal loss, the background under the 
X(3872) — )• TT^ TT~ J / ip signal peak is relatively flat and 
similar in shape to that under the ip' — > TT~^Tr~ J/tp peak. 



MONTE CARLO RESULTS 

We use Monte Carlo (MC) simulated events to deter- 
mine acceptance and to evaluate possible differences in 
mass biases for the ip' and AT (3872) mass regions [3l| . 
The 'ip' MC simulation uses an input mass and width of: 
m^> = 3686.09 MeV and F^. = 0.3 MeV [7]. The defauh 
X(3872) simulation assumes J^^ = 1++ and a Tr'^ir^J/ip 
final state that is entirely B — >■ pJ/ip with the p and 
J/^p in a relative S'-wave [3^- The AT (3872) mesons are 
generated with a mass of M^Z^j^s — 3871.40 MeV and 
zero natural width. The simulated events are processed 
through the same reconstruction and selection codes that 
are used for the real data. 

We perform an unbinned three-dimensional likelihood 
fit (Mbc vs. M{-K+Tr-J/'tp) vs. AE) to the selected 
data using a single Gaussian function for the Mbc signal 
probability density function (PDF) and an ARGUS func- 
tion |33l as the PDF for the combinatorial background 



(i.e., backgrounds where one or more of the tracks used to 
reconstruct the B originates from the accompanying B). 
For AE we use a bifurcated Gaussian for the signal PDF 
and a second-order polynomial for the AE combinatorial 
background. For the M {■k'^ 'k~ J / ip) signal PDF we use a 
Breit-Wigner function (BW) convolved with a resolution 
function that is the sum of a core and tail Gaussian; for 
the combinatorial background PDF we use a third-order 
polynomial. For ip' fits in both data and MC, we fix the 
BW width at 0.3 MeV. For the X(3872) MC fits, we fix 
the BW width at zero. 

In addition to combinatorial background, these crite- 
ria select events of the type B — >■ KxJ /"^P, where Kx 
designates strange meson systems that decay to ATtt+tt" 
final states such as the A'i(1270), A:*(1430), etc. Q. 
The Mbc and AE distributions for these events are the 
same as those of the A" (3872) signal, but they produce a 
slowly varying M{tt'^tt~J/i/j) distribution in the ip' and 
X(3872) signal regions. The Mbc and A^ PDFs that are 
used to represent this peaking background are the same 
as those used for the signal and a linear form is used for 
its M{TT+Tr-J/ip) PDF. 

The results of fits to MC samples of B+ -^ K+ip', 
B" -^ Ksip', B+ ^ K+X{3872) and B" -^ KsX{3872) 
are summarized in Table |ll In order to facilitate com- 
parisons of the resolution for different decay channels, 
the relative fraction of the tail and core Gaussian for all 
modes is fixed at the value returned from the fit to the 
K+ip' MC sample (17.7%). This restriction is found to 
induce negligible differences from the shapes of the res- 
olution functions that are individually optimized for the 
other samples. While the core resolution width is nearly 
the same for all channels, the tail resolution widths for 
A'(3872) decays are significantly higher than those for 
the Ip' , but in both cases the tail widths for the Ar+ and 
Ks modes are consistent with being the same. The MC 
indicates that there are biases in the Af (tt+tt" J/-0) mea- 
surement that are smaller for the AT (3872) modes than for 
the Ip' modes. These are due to a bias in the measure- 
ment of the low momentum charged pions. The pions 
from X(3872) decays have, on average, higher momen- 
tum than those from ip' decays and the X(3872) mass 
measurement bias is smaller. In both cases, the mass 
measurement biases for the AT^ and Ks modes are consis- 
tent with being the same. The results of fits to the com- 
bined A'+ and Ks modes are also shown in Table H] (The 
listed efficiencies do not include the ip' — > 7r+7r~J/'0, 
J/i/j — > i^£^ or Ks — > TT+TT^ branching fractions). 



FITS TO THE iP' -^ tv+tt~J/iP DATA SAMPLES 

For fits to the ip' data we fix the BW width at 0.3 MeV 
and allow the core and tail widths of the M (tt^ ii^ J / ip) 
resolution function to vary as free parameters. The re- 
suhs of the fits to B+ -)■ K+ip' {B° -)■ Ksip') are the 



TABLE I: Results from fits to the selected MC event samples. 
Here e — Nsig/Ng^n is tlie detection efficiency, (Tcoro and ataii 
are the widths of the core and tail components of the mass 



resolution and Mg, 



Mflt are the MC mass measurement 



biases. All errors are statistical. 



Channel 


(percent) 


(MeV) 


0"tail 

(MeV) 


Mgen - Mflt 

(MeV) 


Ksi>' 
Combined 


17.8 ±0.2 
14.1 ±0.2 


1.83 ±0.02 

1.83 ±0.03 

1.84 ±0.02 


5.66 ±0.14 
6.10 ±0.21 
5.66 ±0.13 


0.74 ± 0.02 
0.74 ± 0.03 
0.72 ± 0.02 


K+X{3872) 
KsX(3872) 
Combined 


19.1 ±0.2 

15.2 ±0.2 


1.93 ±0.04 
1.89 ±0.02 
1.93 ±0.02 


7.69 ±0.17 
7.64 ±0.21 

7.70 ±0.15 


0.60 ± 0.02 
0.64 ± 0.02 
0.60 ± 0.02 



smooth curves in the upper (lower) panels of Fig. [H where 
Mbc, M {■K'^ 'K~ J / tji) and AE distributions for events 
within the signal regions of the other two quantities are 
shown. In each panel, the combinatorial background is 
shown as a (red) dotted line, the combinatorial plus peak- 
ing background is shown as a (green) dashed line and the 
total background plus signal is shown as a (blue) solid 
line. The fit results are summarized in Table |IT1 They 
show a mass bias, i.e., a difference between the fitted 
mass and the PDG world-average value for m^/, that is 
larger than the MC mass bias, indicating that the MC 
simulation of the bias in the pion momentum measure- 
ment is imperfect. 



TABLE IL Results from fits to the ip' event candidates. Here 
A^°^*° denotes the number of signal events returned from the 
fit, (Tcorc and crtaii are the mass resolution parameters, and 
AA'/pDG = MpDG — Mat denotes the mass measurement bi- 
ases. All errors are statistical. 



Channel 


^GVtS 


^corc 

(MeV) 


0"tail 

(MeV) 


AMpDG 

(MeV) 


Combined 


3575 ± 64 
814 ± 30 
4367 ± 72 


2.25 ±0.05 
2.45 ±0.11 
2.28 ± 0.04 


8.4 ±0.5 
13.8 ±1.6 
8.7 ±0.5 


1.12 ±0.05 
1.05 ±0.12 
1.11 ±0.05 



As a test of the validity of the MC acceptance calcu- 
lations, we determine branching fractions for B — > K^ip' 
and Ksip' via the relation 



B{B -^ KiIj') 



/xrcvts 



Nbe'^k/kB. 



■0'->ir+ir J/tl>^J/,l,^ee 



(2) 



where iVf '" ^ 



is the number of signal events for K — A'+ 
and K ^ K°, Ngg = (772 ± 11) x 10^ is the num- 
ber of BB events in the data sample, B^i^^+^- j/^ = 
0.336 ± 0.004 and Bj/^^u = 0.119 ± 0.001 (sum of the 
e+e" and iX^ n~ modes) are PDG world-average branch- 
ing fractions [7|, ex is the efficiency for the correspond- 
ing K channel, fx+ = 1 and /a's — 0.346 l35|. The 



results are: B{B+ -^ K+iP') = (6.51 ± 0.12) x lO"-* 
and B{B° -^ K°iP') = (5.22 ± 0.19) x lO"", where only 
statistical errors are shown. The B'^ branching fraction 
result agrees well with the PDG world-average value of 
(6.46 ± 0.33) X 10"^. The B° result is somewhat lower 
than the PDG value of (6.2 ± 0.5) x lO"'' [7], however, 
the errors quoted on the measurements reported here do 
not include systematic uncertainties. 



X(3872) -^ TT+TT- J/V' MASS, WIDTH AND 
PRODUCT BRANCHING FRACTIONS 



The upper panels in Fig. [5] show the Afbc, 
Af (tt+tt^ J/Vj) and AE distributions for events within the 
signal regions of the other two quantities for the i?+ — > 
K'^X{3872) event candidates together with the results 
of the fit. In these fits, the peak mass and full width 
of the BW function that represents the M{tt^tt~J/iIj) 
signal are free parameters, the width of the core Gaus- 
sian resolution function is fixed at cTcoro ~ 2.28 MeV, 
the measured ip' core resolution value, and the width of 
the tail Gaussian is fixed at ataii = 11-5 MeV, the tail 
width from the -0' data sample fit multiplied by the ra- 
tio of the MC-determined X(3872) and ip' tail widths 
to account for its M {tt~^ tt~ J / tp) dependence. The value 



for r 



X(3872) 



returned from the fit is at its lowest allowed 



value of 0.1 MeV [36]. Other results from the fit are 
summarized in Table IIIII 

The lower panels of Fig. [2] show the Afbc, 
Af (tt+tt" J/i/)) and AE distributions for events in the sig- 
nal regions of the other two quantities for the Ks event 
sample, where an X(3872) — )• n'^ ■k~ J / ip signal is evident. 
The results of a fit that fixes the natural width at zero 
and the resolution widths at the same values used for the 
fit to the iiT+X (3872) channel but with the peak mass 
allowed to vary, are shown as curves in the figure and 
summarized in Table IIIII The statistical significance of 
the Ar(3872) signal yield for the Ks event sample is 6.1cr. 
This is determined from — 21n(£o/'Cmax), where £max is 
the maximum likelihood and Cq is the likelihood for zero 
signal yield with the change in the number of degrees of 
freedom taken into account. The difference in mass for 
the X(3872) state produced in 5+ minus that from B^ 
decays {i.e., AM = M+ - Afo) is 

AA/jf(3872) = (-0.69 ± 0.97 (stat) ± 0.19 (syst)) MeV. 

(3) 
Although many sources of systematic error on the mass 
measurement cancel in the this difference, assumptions 
on the natural width used in the fit and possible differ- 
ences in momentum measurement biases between charged 
and neutral kaons do not cancel. We estimate the er- 
ror associated with the natural width to be 0.14 MeV 
from the change in AAfx(3872) determined from a fit to 
the Ks event sample that uses a natural width fixed at 
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FIG. 1: The Mbc (left), M{-k+tv~ J/ip) (center) and AE (right) distributions for B+ -> K+ip' (top) and B° -> Kstp' (bottom) 
event candidates within the signal regions of the other two quantities. The curves show the results of the fits described in the 
text. 



3 MeV. The difference of the measured tp' masses in the 
B+ -J> K+ip' and B° -^ Ksip' channels is AM^/ = 
(-0.07 ± 0.13) MeV. We use the error on AM0/ as an 
estimate of the systematic error associated with possible 
different charged and neutral kaon measurement biases. 
This result strongly disfavors the prediction of 
Ref. [19]. The BaBar measurement for this quantity is 
(2.7 ±1.6 ±0.4) MeV 0. 



TABLE III: Results from fits to the X(3872) event candidates. 
Here jV"^*^ are the numbers of signal events returned from the 
fit and Affjt is the fitted mass value. All errors are statistical. 



Channel 



K+X(2,9,72) 
K°X{3872) 
Combined 



Ar<= 



151 ±15 
21.0 ±5.7 
173 ± 16 



Mfit (MeV) 



3870.85 ± 0.28 
3871.54 ±0.93 
3870.92 ± 0.27 



Mx(3sr2) determination 



Since the mass difference is consistent with zero 
and the resolution functions for the ii'+V(3872) and 
K'^X{3872) are consistent with being the same, we de- 
termine an V(3872) mass value from the single fit to the 



combined samples. To account for the mass measure- 
ment bias, we correct the fitted mass given in Table IIIII 
by adding a correction 5M = (0.92 ± 0.06) MeV, which 
is the MC-determined X(3872) mass measurement bias 
scaled by the ratio of the measured and MC-determined 
■0' mass biases. The validity of this procedure is tested 
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FIG. 2: The Mbc (left), M{tv+tv' J/ip) (center) and AE (right) distributions for B+ -^ K+X{3872) (top) and B° -^ KsX{3872) 
(bottom) event candidates within the signal regions of the other two quantities. The curves show the results of the fit described 
in the text . 



with MC event samples of narrow resonances with ip' 
and X(3872) (J^*^ = 1++) decay dynamics at different 
mass values ranging from m^' to 3872 MeV. It is found 
for both dynamics that the MC mass bias falls linearly 
with increasing M [tt^ -k^ J / ijj) with slopes (6^^'-^) that 
are very nearly equal: h^^^ = -0.096 ± 0.004 keV/MeV 
and b^x(3872) = -0.097 ±0.004 keV/MeV, indicating that 
using the ip' measurement performed at a mass that is 
186 MeV below Afx(3872) to scale the mass shift near 
3872 MeV is reasonable. 

The offset between the MC-determined V^'-like and 
X(3872)-like mass biases is (0.053 ±0.005) MeV. We use 
this offset, scaled by the tp' data-MC mass bias ratio, 
as the systematic error associated with the decay model. 
The systematic error associated with the MC modeling 
of the low energy pion momentum measurements is de- 



termined by comparing results from different versions of 
the MC simulation to be 0.15 MeV. 
The result is 

Mx(3872) = (3871.84 ± 0.27 (stat) ± 0.19 (syst)) MeV, 

(4) 
where the systematic error is dominated by the error on 
the mass bias correction (0.16 MeV) and uncertainties 
in the decay dynamics used to generate the MC sam- 
ples used to study the mass bias (0.09 MeV). It also in- 
cludes the uncertainties in the J/ip and ip' masses and the 
choice of parameterization used in the three dimensional 
fit. The latter is estimated from the quadratic sum of 
the changes induced by ±lcr variations of the fit parame- 
ters and from the use of different functional forms for the 
PDFs. The systematic error evaluation is summarized in 
Table HVl 
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TABLE IV: Systematic errors on the mass measurement. 



Source 


Systematic error (MeV) 


mj/^ 


0.01 


m^' 


0.04 


Bias correction 


0.16 


3-dim. fit model 


0.03 


MC model dependence 


0.09 


Quadrature sum 


0.19 



Jf(3872) 



upper limit 



The current best limit on the width of the X(3872) is 
the 90% confidence level (CL) upper limit of Txi3S72) < 
2.3 MeV reported in the original discovery paper [l|. This 
is narrower than the Af (tt+tt^ J/^) mass resolution of 
the Belle detector in the mass region of the X(3872), 
(a) ~ 4 MeV. However, the three dimensional fits used 
in the analyses reported here are sensitive to natural 
widths that are narrower than the resolution because of 
the constraints on the area of the M [tt'^ ■k~ J / ip) signal 
peak provided by the A/bc and A£' components. Be- 
cause of these constraints on the area of the peak, the 
measured peak height is sensitive to Yx(d.s,72)- This is 
demonstrated in Fig. |3l which shows the results of fits 
to high statistics MC samples where the X(3872) is gen- 
erated with widths ranging from zero to 2.5 MeV. Al- 
though the measurements have some bias, especially at 
very small widths, the different input widths are clearly 
distinguishable. The curve in Fig. [3] shows the results of 
fit of a parabola to the MC measurements. 




1.5 

r Input (MeV) 



FIG. 3: Fitted vaules for rx(3872) (vertical) versus the MC 
generator input values (horizontal). The curve is the result of 
a fit to a second-order polynomial. 

A fit to the X(3872) — > ti^tt^ J/'ip mass peak in data 
with rx(3872) as a free parameter returns a value that 
is at the lower limit imposed on the fit. To establish an 
upper limit on its value, we made a study of how the fit 
likelihood depends on Fx(3872)- 

In the three-dimensional fit, there are correlations be- 
tween the fitted width, the numbers of signal events (jT-sig) 



and peaking background events (ripeak)- The other pa- 
rameters have negligible correlations with the width. We 
therefore performed a series of fits to the data where we 
fixed rx(3872) at a sequence of values ranging from 0.1 to 
3.0 MeV. In these fits all parameters other than ngig and 
"peak were fixed at their best fit values; risig and ripeak 
were allowed to vary. Figure [4] shows how the fit like- 
lihood changes with rx(3872)- The arrow in the figure 



indicates the width value, Fx 



(3872) 



= 0.95 MeV, below 



which 90% of the integrated area under the points is con- 
tained. This value is below the experimental resolution. 
To check sensitivity to uncertainties in the mass resolu- 
tion width, we repeated the scan using the value of the 
tail resolution width determined from fitting the ^' peak 
without any rescaling. This had negligible effect on the 
width of the likelihood. 
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FIG. 4: Likehhood values from the rx(3872) scan described 
in the text. The region of the plot below the arrow contains 
90% of the total area under the points. 

In order to evaluate whether our measured limit is rea- 
sonable given the size of our data sample, we derived 
width upper limits from similar analyses of 24 statisti- 
cally independent, 170-event MC samples that were gen- 
erated with Fx(3872) = 0. Of these, twelve produced 90% 
CL upper limits that are less than 1 MeV; five returned a 
fit value at the lower limit imposed on the fit. In a set of 
24 MC samples generated with Fx(3872) — 1 MeV, none 
returned a width value at the lower limit of the fit and 
17 produced 90% CL lower limits that exclude zero. 

The '0' width has been precisely measured in e^e~ [37[ 
and pp [13 threshold scans to be 0.304 ± 0.009 MeV 0, 
a value that is well below the resolution of our measure- 
ment. We validated our experimental sensitivity to nar- 
row natural widths by refitting the ip' data sample using 
resolution parameters fixed at the values given in Table HIl 
but with F^' left as a free parameter. The fit result is 
F^' — 0.53 ± 0.11 MeV. An examination of the fit like- 
lihood shows that it is well behaved and excludes a zero 
width value by more than 4cr. The measured value is 
0.23 ± 0.11 MeV above the PDG's world-average value, 
which is consistent with the bias value at F ~ 0.3 MeV 
derived from the fitted curve in Fig. [31 namely 0.25 MeV. 

As an upper limit on the natural width of the X(3872), 
we inflate the 90% CL value determined from the scan 



values shown in Fig. |4]by 0.23 MeV, the measured dif- 
ference between our measurement of F^' and its world- 
average value, to account for a possible measurement 
bias. Since both the simulated and observed biases are 
positive and indicate that our measured limit is biased 
high, this produces a conservative value for the upper 
limit. The result is 



F 



X(3872) 



< 1.2 MeV 90%CL, 



(5) 



which is more restrictive than the previous 90% CL limit 
of 2.3 MeV [i]. 



Product branching fractions 



TABLE V: Systematic errors on the product branching frac- 
tion measurement. 



Source 


K+X{2,S72) 


ii'sX(3872) 


Ks/K+ Ratio 




(percent) 


(percent) 


(percent) 


Nbb 


1.4 


1.4 


- 


Secondary BF 


1.0 


1.0 


- 


MC statistics 


1.0 


1.0 


1.4 


MC model 


2.1 


2.1 


- 


Hadron ID 


3.7 


2.6 


1.1 


Lepton ID 


1.1 


1.1 


- 


Tracking 


1.8 


1.4 


0.4 


3-dim. fit model 


3.0 


5.0 


6.0 


Ks efficiency 


- 


4.5 


4.5 


Quadrature sum 


6.0 


8.1 


7.7 



We determine product branching fractions for B^ — > 
K+X, X -^ TT+TT-J/iP and B° -)■ K°X, X -)■ tt+tt'J/^ 
via the relation 



B{B -)■ KX{3872)) x B{X{3872) -^ tt+tt" J/V') 



Arcvts 
^^ K 



^bb^kIk^j/iP- 



(6) 



where the notation is the same as that used for Eq. [2] 
The results are 

B{B+ -^ K+X{3872)) x B{X{3872) -> tt+tt" J/tA) = 
(8.61 ± 0.82 (stat) ± 0.52 (syst)) x 10"^ (7) 



and 



B{B" -^ K°X{3872)) x B{X{3872) 



.+^- 



JM 



(4.3 ± 1.2 (stat) ± 0.4 (syst)) x 10" 



(8) 



where the systematic error includes uncertainties in the 
MC simulation of the tracking, particle identification for 
the leptons and charged kaon, Ks reconstruction, un- 
certainties in the number of BB meson pairs, choice of 
parameterization used in the three dimensional fit, MC 
statistics, decay model dependence and the error on the 
world-average J/ip — ^ £~^£~ branching fraction, all added 
in quadrature. The computations are summarized in Ta- 
ble |Vl The ratio of the B" and B~^ product branching 
fractions is 



RiX) 



B{B° -^ K^X{3872)) 



B{B+ -^ K+X{3872)) 
0.50 ± 0.14 (stat) ± 0.04 (syst). 



(9) 



where the systematic error evaluation is summarized 
in Table |Vl This value is above the range preferred 
by some molecular models for the X(3872): 0.06 < 
R{X) < 0.29 [39]. The BaBar resuh for this ratio is 
rIx) = 0.41 ± 0.24 ± 0.05 ^2Qj]. 



SEARCH FOR A CHARGED PARTNER OF THE 

X(3872) IN B ^ Kir+n^J/tl; DECAYS 

We search for a charged partner of the X(3872) decay- 
ing into Tr~^TT^J/-ip using the selection criteria described 
above for the tt+tt" J/-;/' analysis, with the exception that 
one of the charged pions is replaced by a tt" . For this we 
require two photons with E^ > 35 MeV that reconstruct 
to a tt" — > 77 with a mass-constrained fit x^ ^ 4.0. In the 
event of multiple 7 entries we choose the candidate with 
the best x^ from the tt'^ mass-constrained fit; for multiple 
charged pions, we chose the candidate that produces the 
lowest value of |Ai?|. 

We perform an unbinned two-dimensional (Mbc vs. 
M (tt^ tt'^ J / ip)) maximum likelihood fit to the selected 
event samples using Gaussian and ARGUS function 
PDFs for the Mbc signal and background, and a Crys- 
tal Ball function [40| and third-order polynomial for the 
M(7r+7r°J/^) signal and background, respectively. For 
the peaking background we use the Mbc signal PDF and 
a linear background shape for the M(7r+7r° J/'0) PDF. 
The Crystal Ball function parameters are fixed at val- 
ues returned from fits to samples of Monte Carlo sim- 
ulated B — )• KX~^, X^ — > p'^J/tp events with mx+ = 
3871.7 MeV and Fx+ = 0. The results of the fit to the 
simulated i?° -^ K^X^ sample are shown in the top 
panels of Fig. [5] 

For the data, we do a series of fits with the X^ mass 
restricted to overlapping 10 MeV mass windows covering 
the range 3850 MeV to 3890 MeV. For the K^ X+ chan- 
nel the largest signal yield is 4.2 ± 7.8 events at a mass 
of 3873 ± 6 MeV. The 90% CL upper limit, correspond- 
ing to the signal yield below which 90% of the area of 
the likelihood function is contained, is 17.3 events. For 
the K^X'^ channel, all mass intervals have a zero signal 
yield and the 90% upper limit derived from the likeli- 
hood function for a peak mass fixed at 3873 MeV is 5.4 
events. Mbc and M {tt'^ tt'^ J / ip) plots for the fit to the 
K~ X^ sample with the highest event yield are shown in 
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the middle panel of Fig. [5] The bottom panels of Fig. [5] 
show the results of the fit to the K'^X^ sample with peak 
mass fixed at 3873 MeV. 




f| 

M(J/v|/ rni) (GeV) 

FIG. 5: The Mtc (left), M{-K+-n°J/^) (right) distributions 
for B -^ KX+{3872), X+ -^ p'^J/ip MC events (top) and 
B° ->■ K-TT+Tv°J/ip (middle) and B+ -^ K°TT+Tv°J/ip (bot- 
tom) event candidates in the data, within the signal region 
of the other quantity. The curves show the results of the fits 
described in the text. 

We determine 90% CL product branching fraction up- 
per limits using the relation 



B{B -^ KX+)xB{X+ -^ p^J/i') < 



Arcvts 



iV, 



(10) 
where -^9o%ul i^ fhe upper limit on the event yield for 
each channel, fK+ = 1.0 and fjio — 0.346 (as in Eq. [2|), 
and ex are the MC acceptances reduced by the system- 
atic error: eK+ = 4.5% and ck" = 2.8%. The systematic 
errors are the same as those listed in Table |V] above, 
with the additional inclusion of a 3% systematic error 
associated with data-MC differences in 7r° detection and 



2.5% for the increase in the upper bounds when the res- 
olution parameters of the M {tt'^ tt'^ J / ip) signal PDF are 
varied by ±10%. The systematic errors are 6% for the 
B° -^ K-X+ and 8% for the B+ -^ K°X+ channels. 
The resulting limits are 

BiB" -^ K-X+) X B{X+ -^ p+J/i!) < 4.2 x 10"^ (11) 



and 



B{B+ -^ K°X+) X B(X+ -^ p'^JM) < 6.1 X 10"^ (12) 

The BaBar limits for the same quantities are B{B''' ^• 
K-X+) X B{X+ -^ p'^Jlii) < 5.4 X 10-s and B{B+ -^ 
K"X+) X BiX+ -^ p+J/ip) < 22 X 10-6 '2^. 



ANGULAR CORRELATION STUDIES 

For subsequent analysis, we define a tighter X(3872) 
signal region that extends ±6 MeV around the 
Af (tt+tt" J/-0) signal peak. For background estimates we 
use ±12 MeV sidebands above and below the signal peak 
centered at 3852 MeV and 3892 MeV. There are in total 
165 events in the signal region; the background content, 
determined from the scaled sidebands, is 34 ± 3 events. 

Angular distributions for the sequential decays B — > 
KX{3872), X(3872) -^ pJ/tP, p -^ tt+tt" and J/^p -^ 
l^d'^ for the 1++ and 2 ^ cases are given by the LHCb 
group in Ref. [41]. Since both the B and K mesons 
are scalar particles, an X(3872) meson produced via ex- 
clusive B — )■ KX decays must have a zero component 
of angular momentum along its momentum direction in 
the B rest frame and, thus, its polarization vector, ex, 
must be along this boost direction. This limits the num- 
ber of independent partial-wave amplitudes needed to 
describe the decay. Moreover, angular momentum and 
parity conservation in X(3872) -^ pJ I'P decay implies 
that for 1++ the p and J jip are in an S- and/or D-wa,ye, 
while for 2 ^ they are in a P- and/or _F-wave. Since 
the X(3872) -^ pJ/^J decay occurs at threshold, only the 
lower partial wave in each case is considered. With this 
constraint, the l"*""*" has only one decay amplitude: L = 
and S = 1, where L the p-J/ip orbital angular momen- 
tum and S their spin state. The 2 ^ hypothesis has two 
independent amplitudes: L = \ with S" = 1 or S* = 2, 
which we denote by Bn and B12, respectively. 

We denote by Ox the angle between the J /-ip and the 
direction opposite to the kaon in the X(3872) restframe. 
In the case of .F^ = 1++, the X(3872) -^ pj/i^ de- 
cay produces a p and J /ip in an S'-wave and, thus, the 
distribution in cos6'js: is expected to be flat. For 2 '", 
the final state is P-wave and the cos Ox distribution is 
oc (1 + 3cos^6'x) for B12 — 0, approximately flat for 
|i?ii| ~ I-B12I, and ex sin^ Ox for Bn = 0. For 1++ decays 
to an S'-wave at threshold, the interaction Lagrangian is 



11 




FIG. 6: 

text. 



Definitions of tiie angles x S'nd 9t as described in the 



values (confidence levels) are 3.82 (0.43), 1.76 (0.78) and 
0.56 (0.97) for |cos0x|, |cosx| and |cos0f|, respectively. 




0:2 0:3 0:4 0.5 03 n~ 
IcosOxl 



C\nt o^ ex ■ (ej/i/" ^ ^p)j where ej/^ and ep polarization vec- 
tors. Thus, the three polarization vectors tend to be mu- 
tually perpendicular. In polarized p — > vr+vr^ decays, the 
pions have a cos^ 6 distribution relative to the ep direc- 
tion, while in polarized J/i) — >■ t^t~ , the decay leptons 



have a sin 9 distribution relative to the 



■-.j/i, 



direction. 



To exploit this, we use a coordinate system suggested 
by Rosner ,42] where the x-axis is the direction opposite 
to the kaon (i.e., the ex direction), the x ~ y plane is 
defined by the kaon and 7r+ and the z axis completes a 
right-handed coordinate system. The angle between the 
7r+ direction and the x-axis is designated as x and the 
angle between the ^^ direction and the z-axis as 9^, as 
shown in Fig. |6l In the limit where the J /^ and p are at 
rest in the X rest frame, the expectation for 1++ has the 
distinctive pattern 



(fN 



d cos 9id cos X 



(X sin 9 1 sin x- 



(13) 



The changes in the values of cos x and cos 9e that occur 
when X and 9i are determined in either the J ftp or p rest- 
frames (instead of the X(3872) frame) are much smaller 
than the bin sizes used in this analysis. 

The CDF results on angular correlations used a three- 
dimensional fit to data divided into twelve bins 151. The 



limited statistics of our sample preclude dividing the data 
into enough bins to make a three-dimensional fit feasible. 
Instead we compare one-dimensional histograms of data 
and MC for different hypotheses. 

The data points in Fig. [7] show the | cosx|, | cos^^l and 
|cos6'x| distribution for X(3872) signal region events. 
The dotted histograms indicate the background deter- 
mined from the events in the scaled M(7r+7r^ J/-0) side- 
bands. The solid histogram is the sum of the background 
(dotted histogram) and simulated MC X(3872) — ^ pJ/ip 
events generated with a 1++ (S'-wave only) hypothesis 
and normalized to the observed signal. (The MC sam- 
ples described in this section were generated using the 
partial wave option of EvtGen [32|.) With no other free 
parameters, we find good matches between l"*""^ expec- 
tations and the data for all three distributions: the x^ 
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FIG. 7: The comparisons described in the text for the J = 
l^"*" hypothesis applied to |cosSx| (top), |cosx| (middle) 
and IcosS^I (bottom). The dashed histograms indicate the 
sideband-determined background levels. 

For J^'^ = 2 ^, in addition to the normalization, 
there are two more free parameters that we take to be the 
ratio |Bii|/|i?i2| and the relative phase between Bu and 
Bi2. A comparison of the measured distributions with 
those for a MC simulated 2 '" state with Bu = finds 
poor matches for all three angular distributions: the x^ 
values (confidence levels) are 14.9 (0.005), 48.8 (< 10"'^) 
and 16.5 (0.002) for | cos6'x|, | cosxl and | cos 6*^1, respec- 
tively. For B12 = 0, there are reasonable matches be- 
tween data and MC for the | cosxl (x^ = 6.04, CL=0.20) 
and |cos6'i| (x^ = 1.92, CL=0.75) distribution, but 
poor agreement in the case of the | cos 9x \ comparison 
(X^ = 16.2, CL=0.003). 

We made similar comparisons with simulated event 
samples for a grid of values for |i3ii|/|i3i2| and its rel- 
ative phase. Figure [5] shows the data - MC comparison 
for the case where {B11/B12) — 1.5e^° ', the value for 
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FIG. 8: The comparisons described in the text for the J^*^ — 
2 '" hypothesis appUed to | cos 9x \ (top), | cos x\ (middle) and 
|cos6lf| (bottom) for B11/B12 = 1.5e^''°'. 



which we found the best match. In this case all three 
MC distributions have acceptable x^ values (confidence 
levels): 4.72 (0.32) for cos6lx, 4.60 (0.33) for cosx, and 
5.24 (0.26) for cos0£. The LHCb analysis uses the pa- 



rameter a — 



Bii 



Bii+Bi 



4l|; the values of |i3ii|/|i?i2| and 



the relative phase that are listed above translate into 
a = 0.69e23°*. 

We conclude that with the current level of statistical 
precision we cannot distinguish definitively between the 
1++ and 2 '" assignments. However, while the 2 '" MC 
distributions for all three angles are similar to those for 
1++, they differ in detail, suggesting that in future ex- 
periments with larger data samples, such as LHCb [43J, 
Belle II [4J] and SuperB [45| , three-dimensional fits based 
on the angles discussed here will be able to distinguish 
between the two J^^ hypotheses. 



FITS TO THE M(7r+7r-) DISTRIBUTION 

For even-parity C = -1-1 states the tt+tt" J/ip final state 
would be a p and J/^ primarily in a relative S'-wave, 
while for 2 ^, the p and J/V' would be in a relative P- 
wave. For the S'-wave case, the M(7r"'"7r~) mass distri- 
bution near the upper kinematic limit is modulated by 
the available phase space, which is proportional to k* , the 
J/t/j momentum in the X(3872) restframc. For a J/tp and 
p in a P-wave, the upper boundary is suppressed by an 
additional (fc*)^ centrifugal barrier. Thus, the high-mass 
part of the tt+tt^ invariant mass distribution provides 
some J^ information. 

We extract a background-subtracted M{n^n^) spec- 
trum from a series of two-dimensional (Mbc vs. AE) 
likelihood fits to data in 20 Me V- wide Af(7r+7r^) bins 
covering the range 0.4 GeV < M(7r+7r-) < 0.78 GeV. 
The extracted yields are corrected for the M(7r+7r~)- 
dependence of the experimental acceptance using results 
from simulated data samples of J5 — >• KX, X — !> pJ/V' 
events where the p mass is artificially set at various 
masses with a narrow width. The peaking background re- 
maining in the data is estimated from the M {-k^ 'k~ J / -ip) 
sidebands to be 12 ± 5 events with an M{'K'^Tr~) distri- 
bution that is similar to that of the X(3872) signal. The 
resulting distribution is shown as data points with error 
bars in Fig. [5] 
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FIG. 9: The data points show the background-subtracted, 
relative-efiiciency-corrected M{-k'^tx~) distribution for 
X(3872) — >■ ■k'^ -k' J / ip events. The curves show the results 
of fits using an S-Wave (dashed) and a P-Wave (sohd) BW 
function as described in the text. 



We fit the M{'k^tt ) distribution for events in the 
X(3872) signal region using the parameterization of 
Ref. [l3| 

dN/dm^^ (X {k*f^+'ffx{k*)\BWpim^^)f, (14) 

where k* is defined above, £ is the orbital angular mo- 
mentum value, fox = l.Oa,ndfixik*) = (l-|-P^A:*2)-i/2 
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are Blatt-Weisskopf "barrier factors" 
relativistic BW expression 



'461 and BWp is the 



BWpim, 



V- 



impTp 



(15) 



Here Tp = ^n[g*/go]^["^p/"^7^,^][/lp(g*)//lp(go)]^ where 
g*(TO^^) is the pion momentum in the p rest frame, Qq = 
q*{mp), hp{q)) = (1 + Ry)-^/\. To - 146.2 MeV and 
TOp = 775.5 MeV [;^. The "radii" Rx and Rp are poorly 
known. Generally Rp = 1.5 GeV~^ is used and CDF 
uses values for Rx that are as large as Rx = 5.0 GeV~^. 
(Higher values of Rx reduce the effects of the /c*(2^+i) 
factor and, therefore, make the 5"- and P-wave differences 
smaller.) We take these values as our default settings. 

The smooth curves in Fig. [9] show the results of the 
S'-wave (dashed line) and P-wave (solid line) fits. The S*- 
wave {£ = 0) case fits the data well: x^ /d.o.f. = 17.5/18 
(CL=49%). The P-wave {i = 1) fit is poorer, x^/d.o.f. = 
32.1/18 (CL=2%). Reducing the Blatt-Weisskopf radius 
for the X(3872) makes the P-wave fit worse, increasing 
Rx to 7.0 GeV~^ improves the P-wave fit y^ /d.o.f. to 
26.5/18, which corresponds to a 9.0% CL. Large changes 
in Rp are found to have little effect on the fit quality for 
either case. 

However, both Belle ^] and BaBar [18| have reported 
evidence for the sub-threshold decay process X(3872) — >• 
a; J/f/). The CDF group pointed out that interference be- 
tween the pJ /'>p and ujJ /^ final states, where lo — ^ 7r~'"7r~, 
can have important effect on the M{tt^-k~) lineshape 
near the upper kinematic limit [1^. We therefore re- 
peated the above-described fits with the inclusion of pos- 
sible effects from p-u interference. 

For these fits we use the form given in Eq. [14] with 
PWp(m^,r) replaced by 



BWp^^ cc BWp 



+ r^e'^-BW^ 



(16) 



where BW^ is the same form as BWp with uj meson mass 
and width values substituted for those of the p, r^ is the 
strength of the ut amplitude relative to that of the p, 
and Ajj is their relative phase, which is expected to be 
95° [H. 

We performed fits to the ^^(Tr+Tr") distribution us- 
ing this form weighted by the acceptance with (j),^ fixed 
at 95° and r^ left as a free parameter. Figure [TO] 
shows the results of the S'-wave (dashed line) and P- 
wave (solid line) fits. The inclusion of a small w am- 
plitude {ruj = 0.07 ± 0.05) improves the S'-wave fit to 
X^/d.o.f. = 15.8/17 (54% CL). The P-wave fit returns 
a larger uj contribution, r^ = 0.48lg;^4, and a good fit 
quality: x^ — 14-6 for 17 degrees of freedom {d.o.f.) (62% 
CL). 

The fits have three components: direct p — > tt^tt^ 
(oc |PWpp) and lo — > tt+tt^ (oc r^jPWt^p) contributions 
and a p-oj interference term. The contributions from each 
component for each fit are listed in Table IVIl 
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FIG. 10: The background-subtracted, relative-efficiency- 
corrected M(7r+7r~) distribution for X(3872) — ;■ -K'^-K~J/ip 
events. The curves show the results of fits using an S-Wave 
(dashed line) and a P-Wave (solid line) BW function with 
effects of p-uj interference included. 



TABLE VI: Summary of the results from the p-uj interference 
fit. 



S'-wave 
P-wave 



Nsig 



159 ± 15 
158 ± 15 



0.07 ±0.05 




j+0.20 
'-0.14 



No 



140.9 
93.2 



V„ 



0.6 ±0.5 



No 



17.8 
60.0 



If the low-mass tails of the ui 
TT+TT^ lineshapes are the same [4£ 



— >■ TT^TT tt'^ and oj 
I, we expect 



N{lj -^ tttt) _ B{X{3872) -^ ujj/^) 



■xB{uj 



iVsig 6(V(3872) -^ 77+71- J /ijj) 

(17) 
where the combined result from Belle I47| and BaBar [18[ 
(measured using oj -^ tt+tt^tt^ decays) is B{X{3872) — >• 
LO J /i;)/B{X (3872) ^ 7t+ 7t- J / %jj) = 0.8±0.3. Using this, 
iVsig = 159±15andS(a;^7r+7r-) = 0.0153±0.0013 [7], 
we find an expected value N{uj -^ tt+tt^) ~ 2.0 ± 0.8 
events, which is between the values derived from both 
the S'-wave and P-wave fits and reasonably consistent 
with either case. 



SUMMARY 

We report a measurement of the difference in masses 
of V(3872) mesons produced in P+ -^ K+ 7i+ 71^ J / ijj and 
B° -^ K°7r+7r-J/^ decays, 

AAfx(3872) = (-0.69 ± 0.97 (stat) ± 0.19 (syst)) MeV, 

(18) 
that is consistent with zero and disagrees with theoreti- 
cal predictions based on a diquark-diantiquark model for 
the X(3872) jl^j. We conclude from this that the same 
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particle is produced in the two processes and use a fit to 
the combined neutral and charged B meson data samples 
to determine: 

A^jf(3872) = (3871.84 ±0.27 (stat) ±0.19 (syst)) MeV. 

(19) 
This result agrees with the current PDG world-average 
value of 3871.56±0.22MeV 3 and supersedes Belle's ear- 
lier mass measurement p| , which was based on a 140 fb"^ 
subset of the current data sample. The width of the 
X(3872) signal peak is consistent with the experimental 
mass resolution and we set a 90% CL limit on its natural 
width of rx(3872) < 1-2 MeV, improving on the previous 
limit of 2.3 MeV. 

We report a new measurement of the product branch- 
ing fraction 

B{B+ -^ K+X{3872)) x B{X{3872) -^ tt+tt" J/V-) = 
(8.61 ± 0.82 (stat) ± 0.52 (syst)) x 10"^ (20) 

which supersedes the previous Belle result pj. The 21. 0± 
5.7 signal event yield for B° -^ K°X{3872) translates to 
a ratio of branching fractions 



BiB° ^ K°X{3872)) 
B{B+ ^ K+X{3872)) 



0.50 ±0.14 (stat) ±0.04 (syst). 

(21) 

An examination of the isospin-related B — > 
Kn^TT^J/ip channel shows no evidence for a charged 
partner to the X(3872) decaying as X+ — >■ p^J/ip and we 
determine 90% CL upper limits on the product branching 
fractions B{B -^ KX+) x B(X+ -^ p+J/ip) of 4.2 x 10"^ 
and 6.1 x 10~® for K = K^ and K ~ K'^, respectively, for 
an X^ partner state with mass between 3850 MeV and 
3890 MeV. These limits are well below expectations for 
the X(3872) if it is purely a neutral member of an / = 1 
triplet, in which case decays to the J3 = ±1 partners are 
favored by a factor of two. 

A comparison of angular correlations among the final 
state decay products finds a good match between data 
and MC expectations for J^'^ = 1++ with no free pa- 
rameters (other than the overall normalization). The 
J^'-' = 2 ^ hypothesis has one complex free parame- 
ter and we found a value for which this hypothesis also 
matches the data reasonably well. For this parameter 
value, the differences between 1++ and 2 ^ expectations 
are small but non-zero and a three-dimensional analysis 
based on the angles that we use could distinguish between 
the two cases with the much larger data sets expected at 
the LHCb [ilj. Belle II jiil and SuperB 45] experiments. 
Fits to the M{tt^tt~) mass distribution that only con- 
sider contributions from p — )• tt+tt^ decays favor S'-wave 
(J^ = 1^) over P-wave (J^ = 2~). However, the addi- 
tion of an interfering contribution from isospin- violating 
w — ?> vr+vr^ decays results in acceptable fits for both the 
S'-wave and the P-wave hypotheses. The P-wave fit re- 
quires a more substantial contribution from lo — > tt+tt". 



but with the current limited statistics for X(3872) ^■ 
■K^Ti^J/ip decays and the poor precision on the ratio 
B{X{3872) ^ uj J /il;)/B{X {3872) -^ tt+tt-J/V-), the 
measured uj — > tt^tt" amplitudes that result from fits 
to M(7r+7r~) cannot be used to distinguish between the 
two possibilities. This also may be possible in future ex- 
periments. 
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